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ABSTRACT 

The contribution of the Unresolved Extragalactic Radio Sources to the diffuse 
brightness of the sky was evaluated using the source number - flux measurements 
available in literature. We first optimized the fitting function of the data based 
on number counts distribution. We then computed the brightness temperature 
at various frequencies from 151 MHz to 8440 MHz and derived its spectral de- 
pendence. As expected the frequency dependence can be described by a power 
law with a spectral index 7 ~ —2.7, in agreement with the flux emitted by the 
steep spectrum sources. The contribution of flat spectrum sources becomes rel- 
evant at frequencies above several GHz. Using the data available in literature 
we improved our knowledge of the brightness of the unresolved extragalactic ra- 
dio sources. The results obtained have general validity and they can be used to 
disentangle the various contributions of the sky brightness and to evaluate the 
CMB temperature. 

Subject headings: cosmology: diffuse radiation, radio continuum: galaxies, galax- 
ies: statistics 



1. Introduction 

In the past years measurements of the absolute sky emission at frequency u ^ 1 GHz 
have been carried out to evaluate the brightness temperature {Tcmb) of the Cosmic Microwave 
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Background (CMB). Besides the non-trivial problem of assuring an accurate absolute cali- 
bration of the measured signal, we need to remember that the sky emission is a superposition 
of different contributions. After subtracting the local emissions (mainly due to atmosphere 
inside the main beam, ground and radio frequency interferences in the far side-lobes) the 
sky brightness temperature (Tsky) can be written as: 



T,fcy(z/, a, 6) = Tcmb{i^) + Tgai{u, a, 5) + T; 



UERsiy) 



where Tgai is the emission of our galaxy and Tuers{^) the temperature of the unresolved 
extragalactic radio sources (ITERS). In the present paper we evaluate the UERS brightness 
temperature and its frequency dependence. 

This paper follows a series of others describing the measurements of the sky brightness 
temperature at frequencies close to 1 GHz gathered by th e TRIS experiment wi th an angular 
resolution FWHMtris ~ 20 deg f JZannoni et^bnO?! ). jOervasi et ahlbnOTl ). JXartari et al 



20071 )). The results obtained in the present paper were used to disentangle the components of 
the sky brig htness and to evalua te the CMB temperature at the frequencies z/ = 600, 820 and 
2500 MHz ( flGervasi et al.l 120071 )). The aim of this work is to provide a new estimate of the 
integrated contribution of UERS to the diffuse brightness of the sky. An accurate estimate of 
TuERsi^) is necessary for the TRIS experiment, but also for all the experiments aimed at the 
study of the spectral distortions in the Rayleigh- Jeans tail of the CMB spectrum. Deviations 
from the black-body distribution can be present at low frequency, but the amplitude of the 
distortions at frequencie s around 1 GHz is nowadays constrained by past experiments at the 
level of few tens of mK (IFixen et al.lll996l ). 



Experiments like TRIS (jZannoni et al.l 120071 ) can reach a control of systematics at the 
level of ~50 mK, a remarkable improvement if compared to previous measurements at the 
same frequencies. On the other h and, relying o n the current knowledge of both amplitude 
and spectrum of the UERS signal (lLongairlll966l ). we can estimate that at 600, 820, 1400 and 
2500 MHz (where CMB observations have been carried out in the past) the extra-galactic 
contributio n is respectiv e ly 81 ± 180 mK, 340 ± 80 mK, 79 ± 19 mK and 16 ± 4 niK (se e 
for example ISironi et al.l (jl990l )). Using the current 178 MHz normalization (jLongaiiill966l ). 
for state-of-the-art experiments, this means that the uncertainty associated with the UERS 
at the lowest frequencies (which are the most interesting when looking for CMB spectral 
distortions), is potentially higher than instrumental systematics. In this paper we show that 
by exploiting all the data available in literature we can significantly improve the present 
status of our knowledge about the UERS contribution, and that TRIS-like experiments are 
essentially limited by the current technology. New and updated estimates of the brightness 
temperature of UERS will be useful also for feasibility studies of future experiments in this 
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field. 

This paper is organized as follows: Section[2]discusses the general properties of the UERS 
and the data in literature. In Section [3] we describe the procedure to fit the available number 
counts; in Section H] we calculate the UERS sky brightness and its frequency dependence. 
Finally in Section [5] we discuss the implications of the results obtained for astrophysics and 
cosmology. 



2. The extragalactic radio sources 

2.1. The population of sources 

The unresolved extragalactic radio sources contribute as a blend of point sources to the 
measurements of diffuse emission, especially with poor angular resolution. Actually UERS 
are an inhomogeneous collection of quasars, radio galaxies, and other objects. These can be 
both compact and extended sources with different local radio l uminosity func tion, lifetimes 
and cosmic evolution. An extensive discussion can be found in iLongain (119781 ). 



Usually we can distinguish two populations of radio sources: steep spectrum sources if 
a > 0.5, and flat spectrum sources if a < 0.5, where a is the spectral index of the source 
spectrum (S^u) oc z/~"). Compact radio sources, like quasars, have mostly a flat spectrum 
[a ~ 0) and are most commonly detected at higher frequencies. On the other hand, extended 
sources like radio-galaxies, have a s teep s pectrum (a ~ 0.7 — 0.8) and dominate low frequency 



counts (see ( iPeacock fc Gulllll98ll ) and (ILongairl 119781 )). 



spt 



Steep spectrum sources and flat spectrum sources contribute in different ways to the 
number counts. Flat spectrum sources are important only at high frequency (z/ > 2 GHz). 
In the same range, flat spectrum source counts seem to be comparable to steep spectrum 
source counts for high flu xes, but at low flux es steep spectrum sources a re still dominating. 



as shown for example by ICondonI (119841 ) and iKellermann fc Wall! (119871 ) T he total number 



(1981 


) and 


Danese et al. 


(1987) 



2.2. Isotropy 



The large scale isotropy of the extragalactic radio sources has been studied by [Webster 

(119771 ). He analyzed several samples of sources measured in a cube of 1 Gpc-side, getting an 
upper limit AN/N < 3% on the fluctuation of the source number. This limit is set by the 
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finite statistics of the sources contained in the survey: ~ 10^ 



Franceschini et aLl (119891 ) have evaluated the fluctuation of the source counts assuming 
a poissonian distribution: at = 5 GHz they found fluctuations of the antenna temperature 
AT4/TA < 10"'^ over an angular scale 6* ~ 5 deg. This fluctuation rapidly decreases at larger 
scales. At lower frequencies the fluctuations increase, but we have at most AT^/T^ < 10~^ 
for V = 408 MHz and ATa/Ta < 10"^ for u = 2.5 GHz at an angular scale ^ ~ 20 deg. 

Due to these considerations, the contribution of ITERS to the sky brightness at large 
angular scale is assumed isotropic in the following discussion. Moreover the radio sources are 
supposed to be randomly distributed and the fluctuation in the number counts is assumed 
to be poissonian. A possible anisotropic contribution of UERS is in most case s negligible 
and limited to the region of the super-galactic plane (see (IShaver fc Pierrelll989l )). 



2.3. The data set 



Many source-number vs flux distributions have been produced in the past years: ra- 
dio source counts at low frequency have been published even since the Sixties, while deep 

surveys at higher frequencies have been perform ed recent ly. Con apilation of source count s 

can be found in sever al papers (see for examp le (ICondorull984l ). ( iFranceschini et al.lll989l ). 
( jToffolatti et al.lll998l ). (IWindhorst et al.lll993l )). Most of the data we used can be extracted 
from the references included in these papers. 

We used the counts distributions at the frequencies between 150 and 8000 MHz, spanning 
a frequency range larger than the one covered by the TRIS experiment. In literature we found 
data at eight frequencies: u = 151 MHz, 178 MHz, 408 MHz, 610 MHz, 1.4 GHz, 2.7 GHz, 5.0 
GHz, 8.44 GHz. The complete list of the papers we examined is reported in Table [TJ Some of 
those measurements were performed at a frequency slightly different from the nominal one. 
In this case, the original measurements were scaled to the nominal frequency by assuming a 
dependence of the source flux ^(z/) ~ z/~° ''. Usually the correction is negligible. The number 
counts extracted from Table [1] is shown in Figure [H 



3. Number counts distribution fit 



The fit was performed on the differential number counts normalized to the Euclidean 
distribution of sources E{S) = S^^'^{dN/dS), using the following analytical expression: 
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Q(S) = Q.(S) H- Q.(S) = ; + ; (2) 

The best values of the fit parameters are summarized in Tables |2] and [31 This analytical 
fit function is empirical. It reproduces the distribution of the experimental data with some 
advantages: 1) it is a simple analytical function; 2) it allows to extend the extrapolation 
beyond the available data, because both the amplitude and slope of the tails are well defined, 
both at low and high fluxes; 3) it is built as the sum of two "populations" of sources with 
different emission, but similar behavior and shape; 4) the fitting procedure is independent on 
the source type and evolution but it is applied to the source counts as currently observed. In 
addition, this analytical form is simply a power law, with various indices at different values 
of flux, exactly as the observations suggest. 



3.1. High flux distribution fit 



To get the fit we first considered the source distribution at high fluxes (i.e. we evaluated 
the parameters of the component Qi{S)). For this high flux distribution we have data at 
all the considered frequencies. The best fit parameters are listed in Table [21 We found that 
the values of the spectral indices ei and /?i obtained at all the frequencies are very similar. 
We then decided to take a unique weighted average value for ei and Pi. The parameter 
is particularly well const rained by the data at 1400 MHz, where a large set of data with low 
scatter is available (see (IWhite et al.lll997l ) and Figure [T|). Conversely available data at 178 
MHz do not span a flux range wide enough to constrain the two slopes. 

Source counts at 2700 and 8440 MHz are the less accurate among the full data set 
we used. At both frequencies there are sets of data not completely overlapping and the 
statistics is poor. Therefore, to take into account these uncertainties, we assumed a uniform 
distribution of the parameter Ai. At u = 2700 MHz we took Ai = {6 — 12) x 10~^, while at 
u = 8440 MHz we took Ai = (16-34) x 10""^ (see Table[2]). In this fitting procedure we used 
all the data collected from the papers listed in t he Table [H We exc luded from the fit only 
the seven data points at lowest flux published by [White et al.l (119971 ) at 1400 MHz because 
these data present a roll-off, exactly in the region where a changing of the slope is expected, 
but in the opposite direction. According to the authors, this roll-off is an indication of the 
incompleteness of the survey at the faint-flux limit. 
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3.2. Low flux distribution fit 



We extended the fitting procedure to the counts at low flux, in order to constrain the 
distribution Q2{S). The two parameters A2 and €2 flxed amplitude and slope of the low- 
flux tail of the distribution, while B2 and /?2 are needed to flt the distribution in the region 
showing a change in the slope. Deep counts are available only at 0.61, 1.4, 5 and 8.44 
GHz, but at 8.44 GHz the number of experimental points and their scatter do not allow to 
flt any parameter. In ad dition we considered the model of the low-flux tail published by 
Franceschini et al.l (119891 ) at 1.4 and 5 GHz. This source evolution model flts the data also 
after the addition of the most recent measurements both at 1.4 and 5 GHz. The low-flux tail 
of the model compared with the experimental data and our flt Q{S) are shown in Figure [21 
This evolution model is able to predict accurately the slope of the low-flux tail and w e used 
it to get the value of 62, which is independent on the frequency (IFranceschinil 120071 ). The 
values of 62, obtained at 1.4 and 5 GHz, are fully compatible with the average value of ei 
previously evaluated. Combining the estimates at 1.4 and 5 GHz we get €2 = — 0.856±0.021. 



Since the model by iFranceschini et al.l (119891 ) is not able to flx the amplitude of the 
number counts, i.e. the parameter A2, we evaluated this parameter by using the experimental 
data at 1.4 and 5 GHz. We get: A2/A1 = 0.24 ± 0.02 at 1.4 GHz; A2/A1 = 0.30 ± 0.04 
at 5 GHz. The ratio A2/A1 is almost independent on the frequency. The small change is 
due to the different contribution of the flat-spectrum sources in the total counts at the two 
frequencies in the low flux region (below 10~^ Jy) and in the high flux region (10 mJy - 1 

Jy). 



Using the model of lFranceschini et al.l (119891 ) and lToffolatti et al.l (119981 ) at 1.4, 5 and 8.4 
GHz, we extrapolated the value of the ratio A2/A1 also to the other frequencies. In order to 
evaluate A2/A1 at lower frequencies we estimated the variation of A2/A1 at 1.4 GHz excluding 
from the counts the flat-spectrum sources. This is the extreme condition, which holds at very 
low frequencies. In fact the other sources contributing to the number counts do not change 
signiflcantly with frequency. We obtain in this situation that 0.23 < A2/A1 < 0.24. The same 
result was obtained starting from data and model at 5 GHz. Therefore for the frequencies 
151, 178, 408 and 610 MHz we take the value obtained at 1.4 GHz, but associating a larger 
error bar, due to the uncertainty in the extrapolation procedure: A2/A1 = 0.24 ± 0.04. 

We then e valuated the contribut ion of the flat spectrum sources in the total counts at 
8.44 GHz (see iToffolatti et al.l (119981 )) in comparison with the counts at 5 GHz. In this 
way we estimate at 8.44 GHz the value A2/A1 = 0.31 ± 0.04. At 2.7 GHz we took a value 
constrained by the results obtained at 1.4 and 5 GHz: A2/A1 = 0.24 — 0.30. 



We flnally estimated the B2 and (32- They are important just to deflne the shape of the 
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distribution in the region showing a change in the slope. At 1.4 GHz data are accurate enough 
to constrain both parameters, but (32 can not be constrained at the other frequencies. Since 
the accuracy of these two parameters is not important for the calculation of the integrated 
brightness temperature, we assumed for them the average value for all the frequencies. 

The summary of the best values of all the parameters of Q2{S) is shown in Table O 
The number counts and the function which fit them are shown in Figure [H In conclusion 
we can note that: 1) Ai and Bi, the two frequency dependent parameters of the fit, take 
different values at each frequency. The same is true for the ratio A2/A1. 2) The power law 
indices (ei. Pi, 62 and 132) are frequency independent and we take a common value at all the 
frequencies. 

4. The UERS contribution to the sky diffuse emission 

4.1. Evaluation of the diffuse emission 

The contribution of the UERS {BjjERsiy)) to the sky brightness is evaluated by inte- 
grating the function S{dN/dS) from the largest flux {Smax) of the measured sources down 
to the lowest fluxes {Smin) corresponding to the faintest sources: 

BuERs{i^)= / -j^{v)-SdS (3) 
The brightness temperature Tuers{'^) is by definition: 



TuERsi^) = BuERs{^)7r-;—^ (4) 

2 kB 

ks being the Boltzmann constant. The values of Tuers at the eight frequencies we considered 
are summarized in Table HI 

From the observations we have Smax ~ 10^ Jy (as measured at 151, 408 and 1400 MHz) 
and Smin ~ 10~^ Jy (as measured in the deepest counts at 5 GHz). While sources at higher 
fluxes if present in the surveyed sky region can be easily measured, the limit at low flux 
is set by the confusion limit or by the observation completeness. In other words there is 
no sharp limit at low flux in the population of the sources. We extended the integration 
down to very faint limits, several orders of magnitude below the faintest detected sources 
{Smin ~ 10^^ Jy). When the integration is extended down to Smin ~ 10^^^ Jy the brightness 
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increase by 3 — 4 % and then the value converges. This increment is comparable with the 
total uncertainty we get on the value of the brightness, as shown in Figure [31 

We extended the integration also to higher values of the flux, in order to test also this 
integration limit. Increasing Smax "well beyond the flux of the strongest sources observed, 
the integral change by less than 0.5% and quickly converges. This is a consequence of the 
very low statistics of sources at the highest fluxes. This is a confirmation that the large scale 
brightness is actually not sensitive to the upper limit of integration. 



4.2. Evaluation of the uncertainty 

The error budget takes into account both the fit uncertainties and the number counts 
fluctuations over the observed sky region. The fit uncertainties were evaluated by means of 
Monte Carlo simulations. For each parameter we considered a gaussian distribution with 
standard deviation as reported in Tables [2] and O Only for the values of Ai at 2.7 and 8.44 
GHz and for A2/A1 at 2.7 GHz we assumed a uniform distribution inside the interval reported 
in Tables |2] and O The error bars of the parameters of the fit (in particular A2/A1 and £2) 
include the uncertainty on the extrapolation at the lowest fluxes for the various frequencies. 
We underline that, as shown in Figure [3], the contribution to the brightness temperature of 
the low-flux tail (below ~ 10^^ Jy) is lower than the overall uncertainty reported in Table HI 

The statistical fluctuations of the sources' number counts have no effect in a large 
part of the distribution because the number of sources is quite large. We concentrated on 
the effect of the fluctuation of the few sources with the highest flux. We considered these 
sources randomly distributed and therefore their fluctuation is Poissonian. We evaluated 
the fluctuation of the brightness in a patch of the sky corresponding to the beam of TRIS: 
^TRis ~ 0.1 sr. The upper limit of the contribution to the temperature uncertainty due 
to the fluctuation in the number of sources is directly measured by the maximum of the 
function 



^2 f^'^Jy^sdS 



( Q \ ■J Smin dS /r\ 



plotted in Figure |4] for the specific case of the 1400 MHz data. For all the frequencies 
this maximum falls around ~ 5 Jy, and its value is from 2 to 6 times smaller than the 
corresponding values reported in Table |H Therefore for every frequency, and over the full flux 
range, the error of the brightness temperature is dominated by the statistical uncertainties 
of the fit parameters. 
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The relative error bar of the brightness temperature is 6 — 7% at 151, 408, 610 and 1400 
MHz increasing up to 9% at 5000 MHz. At 178 MHz the available measurements are few 
and old and the obtained error bar is 13%. At 2700 and 8440 MHz both the quantity and 
the quality of the data do not allow accurate estimates of the parameters of the fit and we 
get an uncertainty of 25 — 30%. 



The integrated contribution of the UERS to the brightness of the sky at the various 
frequencies is shown in Figure The distribution can be fitted by a power law: 



Setting z/q = 610 MHz (chosen because it is close to one of the channels of the TRIS experi- 
ment), we obtain the best fit of Tq and 70 shown in TaJole\5\{FITl). As shown in Figure in 
spite of the large error bars at 2700 and 8440 MHz, scatter of the data points is limited. The 
fit of a single power law, done excluding the data with the largest uncertainty (at z/ = 178, 
2700 and 8440 MHz), gives the values of Tq and 70 shown in Tab\e^{FIT2) . Now the scatter 
of the experimental data is much smaller, as shown by the value of the reduced x^- 

In both cases the results obtained are fully compatible with the slope of the steep 
spectrum sources, which are therefore the main contributors to the source counts. However 
it is interesting to check the contribution of the flat spectrum sources especially at high 
frequencies. We assumed a flat spectrum component with fixed slope 71 = —2.00 and a 
dominant steep spectrum component with slope 70 = —2.70: 



The results of the fit of Tq and Ti are shown in Table O {FIT3). Doing so we can get an 
estimate of the contribution of the flat-spectrum component in the number counts at the 
various frequencies. The value of TuERsi^){^/^oy'™ at the frequencies analyzed is shown in 
Figure [6], together with the power law best fits. Figure [6] shows that the two last fits {FIT2 
and FIT3) are equivalent within the error bars. 



4.3. Frequency dependence 




(6) 




(7) 
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5. Discussion 



5.1. Analytical form of the fit function 



As discussed in Section [31 tlie fit function Q{S) is a power law distribution taking 
different slopes and amplitudes in three different regions of the sources' flux. We prefer to 



deal w ith a power law (like Q(S )) instead of a polynomial fit, as performed by lKatgert et al. 



( 1l988l ) and lHopkins et al. A polynomial fit can be better adapted to the experimental 

points, but its validity range is restricted to the region covered by experimental points. 
Therefore it is not possible to use a polynomial fit to extrapolate the distribution outside 
this range. Conversely our fit function can be extrapolated because amplitude and slope of 
the tails are well defined and take into account the shapes expected from the counts evolution 
models (IFranceschini et al.lll989l ). 



According to lLongairl (119781 ) we should expect a broadened maximum in the distribution 
of the source counts with increasing frequency. This indication seem s to be confirra ed looking 
at the differential counts shown in ( (IKellermann fc Walllll987l ) and (ICondonlll984l )). In spite 



of these expectations we have obtained a good fit of the counts at the various frequencies, 
using the same function with the same slopes, as shown Figure [H Part of this effect is 
probably hidden by the larger scatter of the high frequency differential counts. In any case 
the broadening of the maximum could be marginally observed above 1400 MHz and the 
effect is not relevant for the calculation of the contribution to the sky brightness. 



5.2. Frequency dependence 

The spectral dependence of the brightness temperature Tuers follows the expectations 
at low frequency where the number counts are dominated by the sources with steep spectrum: 
a ~ 0.7. At high frequency the situation is more complex. We could expect a flattening, be- 
cause at these frequencies the number counts of flat- spectrum sources begin to be important. 
It will probably appear at frequencies higher than 1400 MHz. Unfortunately the available 
data are not accurate enough to constrain the fit parameters (in particular Ai) at 2700 and 
8440 MHz. The values obtained at these two frequencies, lower than the values expected 
looking at the other frequencies, could also be an indication of incompleteness of the surveys 
(see Figure E]). Conversely at 5000 MHz data are better and numerous. We obtain a value 
for Tuers fully consistent with the slope of the data at low frequency, where steep -spectrum 
sources are dominating, even if there is a marginal indication of a spectral flattening (see 
Figure [6]). In fact when we fit the data adding the contribution of the flat-spectrum sources 
we get an estimate of this contribution which is Ti/Tq ~ 2% at z/ = 610 MHz and Ti/Tq ~ 9% 



at z/ = 5 GHz (see FIT3 in the Table [5]). 



5.3. Previous estimates of Tjjers 



The values of brightness temperature shown in Table H] have error bars of roughly 7% 
at most frequencies. The uncertainty is a bit larger at 178 MHz and even worse at 2700 
and 8440 MHz, because of the quality of the nuinber c o unts data. So fa r very few esti- 



mate s of Tjjers have been published (see (lLongairlll966l ). (IWalllll990l ) and (IBurigana et al. 



2004J )). and the frequencies covered were limited and sometimes the uncertainty was not 
note d. Our results are in agreement with th e values pre viously estimated by iLongair 
1966[ ). Tr/£;^5(178 MHz) = 23 ± 5 K, and bv IWalll Jl99oh . Tr/ER5(408 MHz) ~ 2.6 K, 
TuERsil-"^ GHz) ~ 0.09 K, and TuERsi'^-^ GHz) ~ 0.02 K, but our error bars are defi- 
nitely smaller. The accuracy of the estimated values of Tjjers is particularly important if 
this contribution is to be subtracted to calculate the value of the CMB temperature at low 
frequency. Table H] suggests that the value of Tjjers needs to be accurately evaluated up to 
a frequency of several GHz, because its value is not negligible. 



6. Conclusions 

We used the source number - flux measurements in literature to evaluate the contribution 
of the Unresolved Extragalactic Radio Sources to the diffuse brightness of the sky. We 
analyzed the count distributions at eight frequencies between 150 and 8000 MHz, spanning 



over the frequency range partially covered by the TRIS experiment (see ( iZannoni et al. 



20071)): y = 151 MHz, 178 MHz, 408 MHz, 610 MHz, 1.4 GHz, 2.7 GHz, 5.0 GHz, 8.44 GHz. 



We optimized the fitting function of the experimental number counts distribution. The 
differential number counts {dN/ dS) at the various frequencies are well described by a multi 
power law empirical distribution Q{S) (see Equation [2]) . The amplitudes {Ai and -Bi) are 
frequency dependent parameters of the fit and have different values at each frequency. Con- 
versely the power law indices {ei, f3i, 82 and (32) have a common value at all the frequencies. 

The contribution of the UERS to the sky brightness was evaluated by integrating the 
function S{dN/dS) from the largest flux {Smax = 10^) of the measured sources down to the 
lowest fluxes {Smin = 10~^^) corresponding to the expected faintest sources. We got the 
brightness temperature with a relative error bar of STjjers/Tuers — 6 — 7% at u = 151, 
408, 610 and 1400 MHz, 6Tuers/Tuers - 9% at z/ = 5000 MHz, 6Tuers/Tuers - 13% at 
u = 178 MHz and STuers/Tuers ^ 25 - 30% at z/ = 2700 and 8440 MHz. 



We finally evaluated the spectral dependence of the point source integrated brightness. 
As expected this dependence can be described using a power law with a spectral index 
7o ~ —2.7, in agreement with the frequency dependence of the flux emitted by the steep- 
spectrum sources. We have also tested the contribution of the flat-spectrum sources, adding 
a second component with the slope 71 = —2.0. The contribution of these sources starts to 
be relevant only at frequencies above several GHz. In fact we estimated a contribution by 
flat-spectrum sources ~ 2% at 610 MHz and ~ 9% at 5 GHz. 

The above results were used to evaluate the CMB temperature at frequencies close 
to 1 GHz froni abso lu te measurements o f the sky temperatur e made by our group (see 
dZannoni etaPboOTh . JGervasi et al.ll2007l ). jTartari et ahlboO?! )). 
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Fig. 1. — Data and fit ((5(5*)) of source counts at tlie various frequencies. Data are extracted 
from papers listed in Tabled! At = 1.4 GHz the two terms Q\ and Q2 of the fit are also 
shown. At V = 2.7 and u = 8.44 GHz the fit profiles are obtained using the distribution 
limits of Ai reported in Table [2J 
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Fig. 2. — Data and fit of source counts at v = lA and at = 5 GHz. Tlie low-flux tail 
of the model by iFranceschini et al.l (119891 ) is also shown (grey region). The thickness is an 
indication of the uncertainty in the normalization with the experimental data. 
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Fig. 3. — The integrated brightness temperature at selected frequencies is shown vs the 
lowest flux value of source counts (Sinf). The upper limit of the integration is fixed at 100 
Jy. The region where experimental data are available is also shown (shadowed region). 
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Fig. 4. — Frequency: 1400 MHz. The maximum of the plotted function C{Sjnin) vs S^in 
sets an upper hmit to the contribution of the source number fluctuation to the brightness 
temperature uncertainty in a beam large 0.1 sr (see the text). 
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Fig. 5. — Frequency spectrum of the brightness temperature of the unresolved extragalactic 
radio sources. The continuous hne is FITl (see the text). 
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Fig. 6. — Brightness temperature of the extragalactic radio sources. The values are multi- 
phed by the power law (u/O.GlGHz)^'™ . The fit curves discussed in the text are also shown: 
FITl (dashes); FIT2 (dots); FITS (continuous hue). 
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Table 1. Reference papers for the extragalactic radio source counts. 



Frequency 


Reference 


151 MHz 


Baldwin et al. fl985). Hales et al. fl988). McGilchrist et al. fl990\ 


178 MHz 


Gower fl966). Loneair f 19661. 


408 MHz 


Condon (1984). CoUa et al. (1973). Mills et al. (1973). Pearson & Kus (1978). 




Wall et al. (1980). Birkinshaw (1986). Grueff (1988). Vallee & Roeer (1989). 


610 MHz 


Condon (1984). Wilson & Vallee (1982). Kateert-Merkeliin et al. (1985). 




Kateert (1979). Bondi et al. (2007). 



1.4 GHz Franceschini et al. (1989) . Kateert et al. (1973) . Katgert fc Spinrad (1974) . 



Katgert (1976) . Fomalont et al. (1974) . Machalski (1978 ). Bondi et al. (2003) . 
Condon et al. (1982) . Condon fc Mitchell (1982) . Coleman et al. (1985) . 
Condon fc Mitchell (1984) . Windhorst et al. (1985) . White et al. (1997) . 
Richards (2000) . Prandoni et al. (2001) . Hopkins et al. (2003) . 



2.7 GHz 


Peacock & Wall (1981). Wall et al. (1981). Wall & Peacock (1985). 
Gruppioni et al. (1997). 


5 GHz 


Franceschini et al. (1989). Partridge et al. (1986). Paulinv-Toth et al. (1980). 
Owen et al. (1983). Davis (1971). Bennett et al. (1983). Ledden et al. (1980). 
Fomalont et al. (1984). Paulinv-Toth et al. (1978). Fomalont et al. (1988). 
Prandoni et al. (2006). Ciliegi et al. (2003). 


8.44 GHz 


Toffolatti et al. (1998). Windhorst et al. (1993). Aizu et al. (1987). 
Seielstad (1983). Fomalont et al. (2002). 



Table 2. Estimated values of the parameters of the fit function Qi{S): ei, Pi, Ai and Bi. 



Frequency Parameter Parameter 



(MHz) 






A(xio-^) 


Si(xlO-4) 


151 


-0.95 ± 0.05 


0.37 ±0.06 


1.65 ±0.02 


1.14 ±0.04 


178 




0.31 ±0.04" 


1.9 ±0.3 


1.75 ±0.07 


408 


-0.84 ±0.03 


0.38 ±0.04 


2.60 ±0.05 


4.6 ±0.1 


610 


-0.91 ±0.07 


0.4 ±0.1 


3.04 ± 0.08 


7.9 ±0.6 


1400 


-0.854 ± 0.007 


0.37 ±0.02 


4.65 ±0.02 


23.2 ±0.8 


2700 


-0.76 ± 0.04 


0.37 ± 0.06 


6-12 


60 ±2 


5000 


-0.83 ± 0.05 


0.29 ± 0.06 


8.4 ± 0.4 


105 ±4 


8440 


-0.78 ±0.05 


0.47 ±0.15 


16-34 


138 ±9 


W.Av. 


-0.854 ±0.007 


0.37 ±0.02 







^The value of (3i at 178 MHz has been evaluated using only few data and for 
this reason it is not included in the calculation of the average value. 
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Table 3. Estimated values of the parameters of the fit function Q2{S): e^-, P2, ^2/^1 and 

B2/B,. 



Frequency 
(MHz) 




Parameter 


Parameter 




£9 d-? 




fi2/fii(xl0^) 


151 






0.24 ±0.04" 




178 






0.24 ±0.04" 




408 






0.24 ±0.04" 




610 






0.24 ± 0.04" 


1.14 ±0.4 


1400 


-OA 


^62 ±0.032" 1.47 ±0.15 


0.24 ±0.02 


1.9 ±0.2 


2700 






0.24-0.30 




5000 


-OA 


^52 ±0.027" 


0.30 ±0.04 


2.3 ± 1 


8440 






0.31 ±0.04* 




W.Av. 


-0. 


856 ±0.021 1.47 ±0.15 




1.8 ±0.2 



^Parameter estimated by using the model by iFranceschini et al.l (119891 ). 



'^Parameter estimated by using the model by iToffolatti et al.l (119981 ). 



Table 4. Brightness temperature of the unresolved extragalactic radio sources. 



Frequency (MHz) 


T^ers (mK) 


AT^ers (mK) 


151 


38600 


2300 


178 


22800 


3000 


408 


2600 


200 


610 


910 


70 


1400 


95 


7 


2700 


11 


+4 
-2 


5000 


3.3 


0.3 


8440 


0.41 


+0.15 
-0.06 
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Table 5. Frequency dependence of the brightness temperature of the unresolved 
extragalactic radio sources (see the text) . 





FITl" 


FIT2'' 


FITS^ 


z/Q (MHz) 


610 


610 


610 


To (niK) 


880 ± 28 


903 ± 30 


876 ± 22 


7o 


-2.707 ± 0.027 


-2.681 ±0.029 


-2.70 


Ti (mK) 






18.9 ±0.2 


7i 






-2.00 


Xr 


1.26 


0.11 


0.05 


Q 


0.026 


0.0095 


0.0031 



^FITl has been obtained using all the experimental data and 
fitting a single power law. 

hpjrp2 -j^Qjg been obtained excluding the experimental data at 
u = 178 MHz, 2.7 GHz, 8.44 GHz and fitting a single power law. 

cpjrpg has been obtained excluding the experimental data at 
u = 178 MHz, 2.7 GHz, 8.44 GHz and fitting steep-spectrum plus 
flat-spectrum sources. 



